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Abstract 

By linearizing the inhomogeneous Burgers equation through the Hopf- 
Cole transformation, we formulate the solution of the initial value problem 
of the corresponding linear heat type equation using the Feynman-Kac path 
integral formalism. For illustration, we present the exact solution for the 
forcing term of the form: F(x,t) = u 2 x + f(t). We also present the initial 
value problem solution for the case with a constant forcing term to compare 
with the known result. 
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I. Introduction 



Nonlinear partial differential equations of evolution type such as the KdV, 
nonlinear Schrodinger, and sine-Gordon equations were known to be exactly 
integrable in the sense that there was a transformation which mapped these 
equations into a pair of linear equation that are solvable by the Inverse Scat- 
tering Method 1 . 

These equations are of homogeneous type so the natural question arises 
whether this nice integrability properties also shared by their inhomogeneous 
counterpart. Several attempts had been made to extend the Inverse Scat- 
tering Method to the inhomogeneous case, although no exact solutions was 
found 2 . 

In this paper, we will be considering this problem, but instead of dealing 
with the above nonlinear evolution equations, we choose the much simpler 
one, namely the inhomogeneous Burgers equation, 

du du d 2 u ,_, . . „ „ . 

m +u lE- v 87> =F{x > t:) ' v> ® (L1) 

where the inhomogeneous or source term F (x, t) is a given function. Hence- 
forth, the x and t will be considered as the space and time coordinates, 
respectively. 

Basically, what motivated us to consider the inhomogeneous Burgers 
equation (1.1) was, beside its simplicity, the well known fact that its homoge- 
neous counterpart is linearizable by the celebrated Hopf-Cole transformation 3 , 

I dip . . 

which enables one to solve exactly the corresponding initial value problem. 

Remarkably enough, the Hopf-Cole transformation (1.2) also transform- 
ing the inhomogeneous Burgers equation (1.1) into a linear heat type equation 4,5 , 



dip 
~di 



' d 2 ip 1 „. ' 



ip (1.3a) 



where 

X 

V (x, t)= J dx'F (x', t) + C (t) (1.36) 

c 

with c being an arbitrary constants, while C (t) an arbitrary function of t. 

The outline of the paper is as follows. In sec. II, we summarize the 
general solution of the initial value problem of the inhomogeneous Burgers 
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equation (1.1), in terms of the Feynman-Kac path integral formalism. Sees. 
Ill and IV will be dealing with the explicit exact solutions for special choices 
of the source term: F (x, t) — Q, a constant, and u 2 x+f (t), respectively. We 
also show that our solutions in Sees. Ill and IV in the limit u — > 0, coincide 
with the result one would obtain by simply applying the Orlowsky-Sobczyk 
transformation 6 . 



II. Initial Value Problem and the Feynman-Kac Path 
Integral 

The Hopf-Cole transformation (1.2) with the linear heat equation (1.3) 
enables one to apply the linear method to solve exactly the initial value 
problem (IVP) of the inhomogeneous Burgers equation (1.1). 

In ref. 6, Orlowsky and Sobczyk considered the pure time dependent 
forcing term, F (x, t) — f (t) , and found that the transformations: 

x' = x-<j>(t), t' = t, v! (x, t') =u(x,t)-W it) (2.1a) 

with 

t t 
W(t) = j drf (r) , (t) = J drW (r) (2.16) 
o o 

transform the corresponding inhomogeneous Burgers equation (1.1) to the 
homogeneous form, 

du' ,du' d 2 u' 

Thus, the corresponding initial value problem is trivially solved by applying 
the transformation (2.1). 

Ablowitz and De Lillo then proceeded to consider the forcing term of 
the form F (x, t) — f (t) 5 (x) and solve the corresponding IVP via a linear 
Voltera integral equation 7 . 

In this section, we will formulate the IVP solution for a general source 
term, F(x,t), via the Feynman-Kac path integral method 8 which proceeds 
as follows. 

Let u (x, 0) = u (x) be the prescribed initial value of the inhomogeneous 
Burgers equation (1.1). Then according to the Hopf-Cole transformation 
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(1.2), the corresponding initial value for the linear equation (1.3) is, 



1 y 

ip(y } 0) = exp -— J dx'u' (x') 
zv o 

The solution of the linear equation (1.3) is well-known, namely, 

oo 

iP(x,t) = J dy K(x,t,y,0)^(y,0) 

— oo 

where the kernel K (x, t, y, 0) satisfies the heat kernel equation, 

dK d 2 K 1 



(2.3) 



(2.4) 



dt dx 2 2v 



+ —V(x,t)K = 



(2.5) 



with the initial condition, K (x, 0, y, 0) = 5 (x — y). 

Instead of solving directly the heat kernel equation (2.5), we adopt the 
Feynman-Kac formulation, for which the solution is given in terms of the 
path integral formula: 



K (x, t, y, 0) = J [Dx] exp [S/2v] 
where S is the related action, i.e., 



(2.6) 



z 

S [x (t)} = J dr 



(2.7) 



Following Feynman 8 , we define the path integral (2.6) as, 

/ 1 \ N / 2 °f 

K (x, t, y, 0) = ^lim f- J J dx± . . . dx N exp [— S (xi, . . . , x N ) /2v] 



e^0 



(2.8) 

where e = t n - t„_i, (n = 1, 2, N), t = 0, t N = Ne = t, x n = x (t n ) , 
y = x , x = xjy, and 



N 



S (x, . . . , y, t) = - J2 

Z n=l 



+ v 



X n -\- X n —\ 



) tn- 



(2.9) 



The choice (2.7) for the action S is verified by substituting eq. (2.9) into 
(2.8) and considering the small t — e expansion, in (2.4), to rederive the heat 
equation (1.3) 8 . 
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III. Solutions For a Constant Inhomogeneous Term 



To illustrate the IVP solution generating method that we presented in Sec. 
II, in this section we choose to solve the inhomogeneous Burgers equation 
(1.1) for the simplest case where the source term F (x,t) is chosen to be a 
constant, i.e., F (x,t) = Q. For this case, the corresponding linear equation 
(1.3) is, 



dip 
~dt 



d 2 1 
v^r^ Qxip 



(3.1) 



dx 2 2v 

where in (1.36) we have chosen the lower integral limit, c = 0, and C (t) = 0. 
Throughout this paper we will keep to choose these values of c and C (t) . 

In the following we apply the Feynman-Kac method to solve the corre- 
sponding initial value problem, while the corresponding steady state solution 
is presented in the Appendix A. 

From eq. (1.36) we find that the corresponding action (2.7) is, 



S[(x,t)} = J dr 
o 



1 ( dx s 

2 [d^ 



+ Qx 



(3.2) 



Since the action (3.2) is quadratic, we may simplify the evaluation of the heat 
kernel path integral (2.6) by expanding the action (3.2) around the classical 
path x c i, i.e., 



x (r) = x d (r) + i] (r) 
where x c i obeys the Euler-Lagrange equation, 

d 2 x d 



dr 2 



= Q 



By a straightforward calculation, we obtain 5 , 

1 \ 1/2 



K(x,t,y,0) = 
where S c i is the classical action 



4irvt 



exp [S d /2v] 



(3.3) 



(3.4) 



(3.5a) 



with 



S cl /2v = i{x,y,tf + C{x,t) (3.56) 

2{x-y)- Qt 2 Qt(6x-Qt 2 ) 
Z (x, V, t) = 7^ , C (x, t) = — (3.5c) 



AJvt 
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For an illustration, let us consider the single hump initial condition 3 , 



u (x) = R5 (x) (3.6) 

where 5 (x) is a delta function and R is a constant. The related initial 
condition for the linear equation (3.1), according to (2.3) is, 



ip (x, 0) = exp 



- Yv Rri{x) 



where 



7] [X 



X 

) = lim J dx5 (x) = {°li^ } 



(3.7a) 



(3.76) 



with e > 3 . Substituting (3.5) and (3.7) into (2.4), then the corresponding 
initial value problem solution of the linear heat equation (3.1) is, 



ip(x,t) = 



where £ = £ (x, 0, t) , and 



7T 



-C(x,t) 



2+(e^-l) er/c(£ ) 



2 

erfc (z) = -j= J d£e~ e 



(3.8a) 



(3.86) 



is the complementary error function 9 . 

Substituting (3.8) into the Hopf-Cole transformation (1.2) gives the re- 
quired solution of the inhomogeneous Burgers equation (1.1) with F (x,t) = 
Q, i.e., 

(to {e~ R/2v ~ l) e- ? » 



u (x, t) — Qt + 



(3.9) 



nt 2 + (e- R / 2v - l)erfc(Q 

which obviously reduces to the corresponding homogeneous solution for Q = 
3 . 

Starting from the corresponding homogeneous solution with the single 
hump initial condition (3.6) 3 , one rederives the solution (3.9) directly by 
applying the Orlowsky-Sobczyk transformation (2.1). 



IV. Solutions for F (x, t) = u 2 x + / (t) 

In this section we give an explicit IVP solution for the forcing term 
F(x,t) = uj 2 x + f (t), for which the action (2.7) is still of quadratic type, 
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i.e., 



S[(x,t)} 



dr 



(4.1) 



where u is a constant, and / (£) a given function of time, t. 

Applying the steepest descent method 8 to compute the heat-Kernel K (x, t, y, 0) 
in (2.6) for a quadratic action, as in Sec. Ill, we obtain the exact kernel: 



K(x,t,y, 0) = Nexp[-Sd/2v] 
where Sd is the classical action (see the Appendix B): 



(4.2) 



[(x 2 + y 2 ) coshut — 2xy] 



+ - 



i 



1 _ 

wsinhurf 



5 = - 

c ' 2 sinh u>t 
t t 

x J dr sinh (cut) f (r) + y J dr sinh uj(t — r) / (r) 
o o 

t r 

x j dr j dr' sinh c<j (t — r) sinh (cur') / (r') / (r) 
oo 



(4.3) 



and the corresponding prefactor which is independent of x and y. 

Thus, the corresponding IVP solution for u (x,t), according to the Hopf- 
Cole transformation (1.2), with ip (x,t) given by (2.4), after substituting the 
heat-Kernel (4.2), is: 



u (x, t) = (2Ax + C) 

oo 

J dy ycxp[~A(y+B/2A) 2 /2v]i:(y,0) 

U> — oo 

sinh u)t °? r on 

J dycxp[^A(y+B/2Af/2v\ 



(4.4a) 



where 



B 



A = -uj cothcut 
2 



x + 



sinh ujt sinh 



i r 

— JdTSinhu,(t-T)f(T) 



c = 



sinh cot 



J dr sinhcjt / (r) 



(4.46) 
(4.4c) 

(4.4d) 



For an illustration, let us consider the shock profile initial condition 1 , 

u (x, 0) = (va) [1 - tanh(ax/2)] (4.5) 
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with a a constant, for which the corresponding initial condition for ip(x,t), 
according to (2.3) is, 

ip(x,0) = 2[1 + e- ax ] (4.6) 

Substituting (4.6) into (4.4), then by a straightforward manipulations, we 
obtain the corresponding initial value problem solution: 

2 {aB+va 2 )/A- 

u (x, t) = (2Ax + C) 



2 A sinh tut 



B + 



X -\- e (aB+va 2 )/A 



(4.7) 



which reduces to the celebrated Taylor shock profile solution 1 in the limit: 
u -> 0, and / (t) -> 0. 

Taking instead the limit uj — > 0, but /(£) 7^ 0, we obtain from (4.7), the 
reduced solution, 



u(x, t) = W(t) + (va) [1 - tanh ({a(x + <p(t) - at} /2)] 



(4.8) 



where <p(t) and W(t) are given by (2.16), as one would expect by applying 
the Orlowsky-Sobczyk transformation (2.1) to the I VP solution of the corre- 
sponding homogeneous Burgers equation 1 ' 6 . It should be mentioned that in 
obtaining the solution (4.8), we have used the following limiting values of A, 
B, and C, obtainable from (4.4), i.e, 



A =2t 



B = -- t +\ t Jdr{t-T)f{r) = -^ + <P{t) 



C = - 

t 



drf{r) 



L0 



where in (4.96), we use the fact that: 



(4.9a) 
(4.96) 

(4.9c) 



C IT I 

cj)(t) = J drW(T) = jdr J dr'f(r') = J dr'(t - r')/(r') (4.10) 
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Appendix A. The steady state solution for F (x, t) = Q 

In the following, we present the steady state solution of the heat equation 
(3.1). First of all, let us look for a single state solution, for which tp (x,t) is 
given by, 

rl>(x,t)=T(t)<p(x) (A.l) 
Substituting (A.l) into (3.1) gives the separated equations: 

<n cuT (A.2) 



dt 

d 2 Lfl 

dx 2 

where 



zuj = (A.3) 



-i+^-®'-'-(tT (a ' 4) 

and u is the corresponding separation constant. 
The solution of eq. (A.3) is well-known, i.e.: 

up (x) = [aAi (z) + bBi (z)} (A.5) 

where Ai (z) and Bi (z) are the Airy functions of the first and second kind, 
respectively 9 , with a and b are constants which are fixed by the boundary 
conditions. 

Thus, the single state solution of eq. (3.1) is given by 

i) (x, t) = [aAi (z) + bBi (z)} (A.6) 

Unfortunately, the corresponding solution for u (x, t) under the Hopf-Cole 
transformation (1.2), is a static solution. 

However, since eq. (3.1) is linear, we may take the superposition, 

N>1 

i> (x, t) = J2 eWnt KAi (z n ) + b n Bi (z n )\ (A.7a) 

n=l 



where 

y + A n , X n = (^J I 2 (A.76) 

as the solution for eq. (3.1). The corresponding solution for the inhomo- 
geneous Burgers equation (1.1), with F (x,t) = Q, through the Hopf-Cole 
transformation (1.2) is, 

N>1 

9v E \a n Ai> (z n ) + b n Bi' {z n )\ 

u (x, t) = --^L (A.8) 

E e" nt KAi (z n ) + b n Bi {z n )\ 

71 = 1 



with Ai' (z n ) and Bi' (z n ) are the derivatives of the corresponding Airy func- 
tions with respect to z n . The solution (A. 8) is obviously non-static. 

Appendix B. Derivation of the classical action (4.3) 

In this Appendix, we present the derivation of the classical action (4.3). 
For this purpose, let us rewrite the action (4.1) in the following form: 

dx t =t f f , j d 2 x 1 2 2 



2S = x (r) — |Sf -Jdrix-^- -wV - 2f (r) x (B.l) 



u 



Then, by using the equation of motion: 

d 2 x 



<h2 u 2 x fir) (B.2) 

the corresponding classical action is, 

, */ 

25 = x(r)^|Sf- Jdrx(r)f(r) (B.3) 

U 

where x (r) is given by the solution of (B.2), with the boundary conditions, 
x (ti) = Xi and x (tf) = Xf. 

Using the Green's function method, the solution of (B.2) is, 

*/ 

x(r) = j dr'G (r', r) f (r) + x (r') -^G (r', r) |Sf (B.4) 

where the Green's function G satisfies: 

^- 2 G (r', r) + u 2 G (r', r) = 5 (r' - r) (B.5) 

with the boundary conditions, 

G(r',t i )=G(r',t f ) = (B.6) 

Solving eq. (B.5) for r < r' and r > r', then using the boundary conditions 
(B.6), the continuity condition, 

G (t' + e, t') —G(t' — e, r') — 0, (e -> 0) (B.7) 
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and the discontinuity condition, 

(r, t') \ t=t>+£ -j-G (r, t') \ t=t ,. £ = 0, (e -> 0) (B.8) 

we obtain: 

q r i\ _ 1 j sinhcj (r' — t/) sinhcj (t — Tj) , r < r' 1 

wsinhcj(r / -r i )\ sinhw(r / -r i )sinhcj(r-r / ), r>r' J 

(B.9) 

Substituting (B.9) in (B.4) gives: 

X W = sinh^-r.) ^ Sillh U ( T ~ ^ ~ ^ Sillh W ( T ~ T /)] 



— w <^ 7 ~'' sinh (t — tj) sinhw (r' — Tj) / (r') 



(B.10) 



— - / dr' sinhcj (r — Tj) sinhw (r' — t/) / (r') 



Using (B.10) in (B.3), and choose :Ti — 0, Tf — t, Xi — y, Xf — x, then by a 
straightforward manipulations we obtain the classical action in (4.3). 
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